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Abstract. Predominantly preliminary single and double inclusive momentum 
spectra of 158 AGeV Pb+Pb collisions, recently measured by the NA44 and 
NA49 Collaborations, are reproduced using the relativistic hydrodynamical 
model HYLANDER-C. Two different equations of state, which both contain 
a phase transition to a quark-gluon plasma, can be used to reproduce the 
(preliminary) data. The space-time geometries in the two calculations differ 
strongly. However, the Bose-Einstein correlation functions of identical pion 
pairs do not show such a strong, but still a significant sensitivity to the effects 
of the equations of state. 


1. Introduction 

Fluid dynamics provides an intuitively simple description of heavy-ion collisions: 
two nuclei smash into each other and are rapidly thermalized and compressed; 
the resulting zone of very hot and dense nuclear matter (the so-called fireball) 
then expands and breaks up into bits of hadronic matter that ultimately reach the 
detectors. In a 158 AGeV Pb-I-Pb reaction the velocity of the incoming projectile 
is 0.99998 c; since the sound speed in the ground state of nuclear matter is about ~ 
0.2 c, shock waves are formed. The density achieved in the center of the hot fireball 
is calculated to be as much as 20 times that of normal nuclear matter. Under 
such extreme conditions, a quark-gluon plasma (QGP) can be formed, which in our 
calculations can reach an appreciable fraction of the total matter [ |^,[ ^. 

Many experimental observables have been suggested with which one might 
ascertain whether a quark-gluon plasma (QGP) has been formed in relativistic 
heavy-ion collisions, or with which one can obtain some idea of the evolution of 
the experimentally generated fireballs. By measuring the correlations of identical 
particles (Bose-Einstein correlations (BEG)), one can get a measure of the sizes 
and lifetimes of the fireballs [ This is very similar to the concept behind the 
Hanbury-Brown/Twiss effect used to measure the size of stars [|^. 
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Fig. 1. Pressure, P, entropy density, s, energy density, e, ratio of pressure and 
energy density. P/e, speed of sound, Cq, and temperature, T, as functions of T 
and/or e, for the equations of state EOS-I (solid lines) and EOS-II (dashed lines), 
respectively (see text). 
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It is the purpose of this paper to use in a hydrodynamical analysis two not too 
different equations of state (EOS) in the attempt to reproduce the predominantly 
preliminary single inclusive momentum spectra of 158 AGeV Pb+Pb collisions, 
recently measured by the NA44 [||],[ H and NA49 Collaborations [0 , and to learn 
which details of the calculations are most sensitive to the effects of the EOS. We 
shall also learn about the initial conditions which one has to use while considering 
the chosen equations of state. After a comparison of the theoretical calculations 
with the measurements, we exhibit the space-time geometries of the real hadron 
sources in the two scenarios, and discuss their particular features. Finally, we shall 
calculate inverse width parameters of Bose-Einstein correlation functions of identical 
pion pairs for the two different hydrodynamical solutions and compare the results 
to preliminary BEC data taken by the NA49 Collaboration [ 

2. The Model, the Equations of State, and the Initial Condi¬ 
tions 

2.1. The Model 

Among the large number of models (cf. refs. [ |.[|] and refs, therein) which 
apply relativistic hydrodynamics to relativistic heavy-ion collisions, we choose for 
the following to use the model^ HYLANDER-C. 

This model applies 3-1-1-dimensional relativistic one-fluid-dynamics, and pro¬ 
vides fully three-dimensional solutions of the hydrodynamical relativistic Euler- 
equations [ 0 - One has to specify an equation of state, initial distributions, e.g., 
in form of parameterizations with initial parameters {cf., refs. [ 0), and a 

break-up (freeze-out) condition in order to obtain an unambiguous solution from 
the hydrodynamical equations. In the calculations we shall assume that hadroniza- 
tion occurs for all particle species at the same fixed freeze-out temperature Tf = 
139 MeV. 

All model calculations are based on the assumption of thermal as well as on 
chemical equilibrium. In both types of spectra (single inclusive momentum spec¬ 
tra and BEC, respectively) we shall include resonance decays. In particular, the 
momentum distributions are calculated in terms of the generalized Cooper-Frye for¬ 
mula (see ref. [0), where explicitly a baryon and a strangeness chemical potential 
are taken into account [ 0. The subsequent calculations of Bose-Einstein correla¬ 
tions are performed using the formalism outlined in ref. [ 0 . The hadron source 
is assumed to be fully chaotic; the influence of partial coherence [ 0 will not be 
considered here. 

2.2. The Equations of State 

Any type of EOS can be considered when solving the relativistic Euler-equations. 
In Fig. I the two equations of state, which we are going to use in the following, are 
plotted in many different representations. 
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Fig. 2. Initial distributions of the energy density, e, the baryon density, ub, and 
the fluid rapidity, y, plotted against the longitudinal coordinate z. The solid lines 
indicate the initial distributions for the use with EOS-I, whereas the dashed lines 
indicate the initial distributions for the use with EOS-II, respectively. 


The first equation of state, EOS-I, has been used in many recent calculations 
(c/., e.g. refs. [ |,[|,[|ni,[[|), in the successful attempt to reproduce hadronic 
single inclusive momentum spectra and BEG, which have been measured in relativis¬ 
tic heavy-ion collision experiments. In particular, EOS-I exhibits a phase transition 
to a quark-gluon plasma at a critical temperature Tc = 200 MeV with a critical 
energy density Cc = 3.2 GeVjfw? (c/., refs. [ @,[[n ,[ |^). 

The second equation of state, EOS-II, is also a lattice QCD based EOS which 
has recently become very popular in the field of relativistic heavy-ion physics (c/., 
ref. [Hi)- This equation of state includes a phase transition to a quark-qluon 
plasma at = 160 MeV with a critical energy density Cc ~ 1-5 GeV/fm?. 

Both EOS have no dependence on the baryon density. For instance, in Fig. 1 
the plot of P{e)/e emphasizes the existence of a minimum P/e at e = Cc = 3.2 
GeV/fm? (w 1.5 GeV/fm?) for EOS-I (EOS-II), referred to as the softest point of 
the EOS. It corresponds to the boundary between the generalized mixed phase and 
the QGP [|^. As it can be seen in Fig. 1, EOS-II yields a much softer equation of 
state than EOS-I. 


2.3. The Initial Conditions 

In the following it is assumed that due to an experimental uncertainty for the 
centrality of the collision, only 90% of the total available energy and the total 
baryon number have been observed. It is then possible to find initial distributions 
for the two here considered equations of state, such that one can reproduce the 
(preliminary) single inclusive momentum spectra of 158 AGeV Pb-fPb collisions, 
which have been measured recently by the NA44 and NA49 Gollaborations. 
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Fig. 2 shows the initial distributions for the energy density, e{z), the baryon 
density, ub^z), and the fluid rapidity, y{z), plotted against the longitudinal coor¬ 
dinate z. We use here for the initial distributions the initial condition scenarios 
which have been extensively described in refs. [ and [ |^. In particular, it is 
assumed that an initial transverse fluid velocity field is completely absent, and the 
initial longitudinal distributions for energy density, e, and baryon density, ub, are 
smeared out with a Woods-Saxon parametrization in the transverse direction, rj_ 
{cf., ref. [ 1^). The choices for the initial parameters, which in each case provide 
the best fit results for the hadronic single inclusive momentum spectra of the two 
considered experiments, are in case of EOS-I (EOS-II): relative fraction of thermal 
energy in the central fireball, = 0.55 (0.20), longitudinal extension of the central 
fireball, A = 1.20 fm (1.00 fm), rapidity at the edge of the central fireball, j/a = 
1.00 (0.85), rapidity at maximum of initial baryon y distribution, y^ = 0.80 (1.50), 
and width of initial baryon y distribution, a — 0.32 (1.00), respectively. 

The maximum initial energy density is ca = 15.3 GeVjfw? (6.55 GeV/ 
and the maximum initial baryon density is = 3.93 frn~^ (1.24 for 

EOS-I (EOS-II), respectively. 71% (30%) of the baryonic matter is initially located 
in the central fireball region. Hence, the use of EOS-I predicts a much larger 
stopping than the use of EOS-II. 

It should be stressed here, that it is really neccessary to use such extremely 
different initial conditions for EOS-I and EOS-II. If one just uses different equations 
of state without changing the initial conditions, which in fact are unknown for 
158 AGeV PbH-Pb collisions, one would have obtained single inclusive momentum 
spectra which differ by up to 20% - 30% when comparing the calculations for EOS-I 
and EOS-II. 

3. Single Inclusive Momentum Spectra 

Figs. 3-6 show the results of the hydrodynamical calculations compared to the 
preliminary and final single inclusive momentum spectra of the 158 AGeV Pb-fPb 
collisions, which have been measured by the NA44 [ ^,[ ^ and NA49 Collaborations 
[ All single inclusive momentum spectra have been evaluated in the nucleus- 
nucleus center of mass system {ycm = 2.91). 

In case of EOS-I (EOS-II) the average value for the baryonic chemical poten¬ 
tial is (/is) = 324 MeV (360 MeV), and the average value for the strangeness 
chemical potential is {fj,s) = 55 MeV (69 MeV), using Tj = 139 MeV for the 
freeze-out temperature. The energy density at freeze-out is e/ = 0.292 GeV/frrA 
(0.126 GeV/frrA). Both calculations yield considerably good agreement with both 
experiments. However, the calculations for the much softer EOS-II already gives 
larger slopes in the transverse mass spectra compared to the calculations for EOS-I. 
The different slopes in the transverse mass spectra have their origin in the different 
transverse velocity fields at freeze-out. For EOS-I we obtain a maximum transverse 
velocity u™“^(/) = 0.46 c, whereas for EOS-II we get the much smaller value 
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Fig. 3. Rapidity spectra, dN/dy, and transverse mass spectra, l/m±dN/dm±, for 
negative hadrons, h~ [i.e., 7r“, K~, p), net protons (without contributions from 
decay), p—p, and neutral kaons, Kg, respectively. The solid lines indicate the results 
of the calculations when using equation of state EOS-I, whereas the dashed lines 
represent the results when using equation of state EOS-II. The shown preliminary 
data have been taken by the NA49 Collaboration [ 
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Fig. 4. Transverse mass spectra, l/m±dN/dm±, for pions, and tt , kaons, 
K'^ and K~, protons (without contributions from A*' decay), p, and anti-protons, 
p, respectively. The solid lines indicate the results of the calculations when using 
equation of state EOS-I, whereas the dashed lines represent the results when using 
equation of state EOS-II. The shown preliminary data have been taken by the NA49 
Collaboration [ ^ . 
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Fig. 5. Unnormalized transverse mass spectra, l/m±dN/dm±, for pions, tt’*' and 
TT~, kaons, and K~, protons (including contributions from A° decay), p, and 
anti-protons, p, respectively. The solid lines indicate the results of the calculations 
when using equation of state EOS-I, whereas the dashed lines represent the results 
when using equation of state EOS-II. The shown preliminary data have been taken 
by the NA44 Collaboration [ . 
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Fig. 6. Rapidity spectra, dN/dy, and transverse mass spectra, l/m±dN/dmj_, 
for protons, p, respectively. The solid (dotted) lines indicate the results of the 
calculations when using equation of state EOS-I including (without) contributions 
from A° decay, whereas the dashed (dashed-dotted) lines represent the results when 
using equation of state EOS-II including (without) contributions from A° decay. 
The data have been taken by the NA44 Collaboration [ 

^max(jj) ^ Q 3Q 

Since the data analysis for the 158 AGeV Pb-|-Pb collisions has not been finished 
yet by the NA44 and NA49 Collaborations, we cannot draw too many conclusions 
here. Therefore, in the following we shall look into the space-time features of the 
hydrodynamical model calculations and try to learn, whether Bose-Einstein corre¬ 
lations provide a more distinct observable in the attempt to decide which one of the 
two considered EOS gives a better description of the data. 

4. Space-Time Geometry and Bose-Einstein Correlations 

In Fig. 7 it can be seen, that the hydrodynamical solutions which are provided from 
the numerical analysis represent an evolution of initally disk-shaped fireballs which 
emit hadrons from the very beginning of their formation. While the relativistic fluids 
expand in longitudinal and in transverse directions, the longitudinal positions of the 
freeze-out points increase their distance relative to the center. Because of the effect 
of transverse inwardly moving rarefaction waves, the transverse freeze-out positions 
move towards the center of the fireball. In the late stage of the hydrodynamic 
expansion the hadron-emitting fireballs separate into two parts while cooling down 
until they cease to emit. 

For both equations of state we obtain fireballs which have similar transverse 
sizes and a QGP phase of similar lifetime. For EOS-I (EOS-II) the total lifetime of 
the QGP is tgcp = 2.4 fm/c (3.0 fm/c), considering the critical temperature Tc 
— 200 MeV (160 MeV). However, the softer EOS-II results in a fireball that has a 
much larger longitudinal size than the fireball which is governed by EOS-I. 
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Fig. 7. Time contour plots of the freeze-out hypersurfaces in the z — r_L plane, 
and temperature evolution plots at r± = 0, for the relativistic fluids governed by 
EOS-I and EOS-II, respectively. In the time contour plots, each line represents 
an isotherme {Tf = 140 MeV) at a fixed time t (timesteps At = 1.0 fm/c). In 
the temperature evolution plots, each line corresponds to a fixed temperature, T 
(temperature steps AT = 20 MeV). 

The total lifetime of the fireball is t^ax = 29.3 fm/c (13.1 fm/c) for EOS-II 
(EOS-I). Such large differences (factor ^ 2) in the lifetimes and longitudinal sizes of 
the two considered scenarios should show up in the inverse widths of the transverse 
“out” [ ^ and longitudinal Bose-Einstein correlations of identical pion pairs. 

In Fig. 8 inverse widths of BEG functions of identical negative pion pairs (in¬ 
cluding decay contributions from resonances) are shown, which have been extracted 
with a Gaussian fit as explained in refs. [ plj,[ |^ , in comparison with preliminary 
measurements of the NA49 Collaboration [ . The inverse widths are plotted as 

functions of the average transverse momentum of the pion pair, K±. The effec¬ 
tive longitudinal radii, Riong, are evaluated in the longitudinal comoving system 
(LCMS), whereas all other effective transverse radii, Rside and Rout, have been cal¬ 
culated in the nucleus-nucleus center of mass system. Since BEG of negative pions 
are in general not of Gaussian shape, (c/., e.g., refs. [ |^,[ |^,[ ^), the fit to a 





















Extracting the Equation of State through Hydrodynamical Analysis 


11 


Gaussian generates an error for the inverse width parameters. 

Therefore, in Fig. 8 also the true inverse widths of the correlation functions 
are included at their 68% level, i.e., the inverse widths at 0.68 • {Iq — 1), where 
lo = C 2 {k, k) is the intercept of the BEG function, (72(^1, fe), of two identical pions 
with momenta fci and ^ 2 , respectively. 

As expected, the transverse radii Rside as functions of A'x are very similar 
when comparing the calculations for EOS-I with the calculations for EOS-II. The 
transverse radii Rout and the longitudinal radii Riong are larger for the calculation 
using EOS-II compared to the calculation using EOS-I. But the differences, e.g., 
in the longitudinal effective radii are not so pronounced anymore as in the case of 
the longitudinal sizes of the exhibited space-time geometries of the hadron sources. 
However, they are quite significant. 

Unfortunately, the data analysis of the Bose-Einstein correlations for the 158 
AGeV Pb-|-Pb collisions has not been finished yet either by the NA49 Collaboration, 
so that the experimental situation does not provide an definite answer yet, which 
one of the considered EOS gives a better description of the data. 

5. Conclusions 

The predominantly preliminary single and double inclusive momentum spectra of 
158 AGeV Pb-fPb collisions, recently measured by the NA44 and NA49 Col¬ 
laborations, have been reproduced using the relativistic hydrodynamical model 
HYLANDER-C. However, the calculations for the much softer EOS-II yield larger 
slopes in many of the hadronic transverse mass spectra compared to the calculations 
for EOS-I. This effect is particularly strong in the transverse mass spectra of protons 
and anti-protons. The initial conditions, which have been found in the numerical 
analysis, differ strongly for the two considered scenarios. The space-time geometries 
in longitudinal direction and in time show a difference in the two calculations of 
more than a factor of 2. In contrast to the common belief [ |^, the Bose-Einstein 
correlation functions of identical pion pairs do not show such a strong, but still a 
significant sensitivity to the effects of the equations of state. 
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Fig. 8. Inverse widths of BEC functions of identical negative pion pairs (including 
decay contributions from resonances) as functions of the average transverse momen¬ 
tum of the boson pair, K±, in the indicated ranges of the particle pair rapidities, 
PK, compared to preliminary data of the NA49 Collaboration [||], respectively. The 
solid (dashed) lines indicate the inverse widths of BEC functions extracted from a 
Gaussian fit for the calculation using EOS-I (EOS-II). The dotted (dashed-dotted) 
lines are the true inverse widths of the correlation functions at their 68% level (see 
text), and the grey zones reflect the theoretical uncertainties, when extracting the 
inverse widths. 
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b. HYLANDER-C is an improved version of the original model HYLANDER {cf., 
refs. 
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